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Abstract. The interaction of mesoscopic interference devices with nonclassical 
electromagnetic fields is studied. The external quantum fields induce a phase factor 
on the electric charges. This phase factor, which is a generalization of the standard 
Aharonov-Bohm phase factor, is in the case of nonclassical electromagnetic fields a 
quantum mechanical operator. Its expectation value depends on the density matrix 
describing the nonclassical photons and determines the interference. Several examples 
are discussed, which show that the quantum noise of the nonclassical photons destroys 
slightly the electron interference fringes. An interesting application arises in the context 
of distant electron interference devices, irradiated with entangled photons. In this case 
the interfering electrons in the two devices become entangled. The same ideas are 
applied in the context of SQUID rings irradiated with nonclassical electromagnetic 
fields. It is shown that the statistics of the Cooper pairs tunneling through the 
Josephson junction depend on the statistics of the photons. 
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1. Introduction 

Interference of electrons in the presence of a magnetostatic flux has been studied for a 
long time since the work of Aharonov and Bohm [T] . The Aharonov-Bohm phase factor 
is acquired by electric charges that encircle a magnetic flux, even if the flux vanishes in 
the vicinity of the paths of the charges. The effect has inspired numerous applications in 
solid state physics 2\. In particular we mention extensive theoretical and experimental 
research on persistent currents in mesoscopic rings [HI El El El E] • 

Electron interference in the presence of a time-dependent magnetic flux (i.e., 
electromagnetic fields) has also been studied [HI EH HU E21 EH Hlj • The intention here is 
not to prove the reality of the vector potential, but to study how electromagnetic fields 
affect interfering electrons. 

The next step in this line of research is to consider nonclassical electromagnetic fields 
[Tallin] which are carefully prepared in a particular quantum state, and study their effect 
on quantum interference [01 CDj . In this case it is shown that the quantum noise in the 
electromagnetic field destroys partly the electron interference fringes. Different types 
of nonclassical electromagnetic fields are characterized by different quantum statistics; 
and we will show explicitly that the electron interference results depend on the photon 
statistics. 

An important feature of multimode quantum electromagnetic fields is entanglement. 
Two electromagnetic field modes can be factorizable (uncorrelated); or separable 
(classically correlated); or entangled (quantum mechanically correlated) [E|. Entangled 
electromagnetic fields have been produced experimentally in laboratories for a long time 
[TH] . In the context of this review article, we consider two distant mesoscopic electron 
interference devices that are irradiated with a two-mode nonclassical electromagnetic 
field [10]. Each field mode is coupled to one of the mesoscopic devices. For entangled 
electromagnetic fields, the electric currents in the distant mesoscopic devices become 
correlated. Moreover the induced correlations of the electrons depend on the nature of 
the correlation between the external photons. 

Similar phenomena can be studied in the context of superconducting quantum 
interference devices (SQUID) [^IIllSJESEil!- Experimental work has so far concentrated 
on the interaction of mesoscopic devices with classical electromagnetic fields. However 
the interaction of a Josephson device with a single microwave photon has recently been 
studied experimentally in reference [21] . 

The interaction of mesoscopic SQUID rings with nonclassical electromagnetic fields 
has been studied theoretically in [23 EE1 EH EE] • In this case the Josephson current is 
a quantum mechanical operator, whose expectation value with respect to the density 
matrix of the external photons, yields the observed current. The interaction of entangled 
electromagnetic fields with two spatially separated SQUID rings has been studied in 
|29| ISO] . It has been shown that the photon correlations can be transferred to the 
Cooper pair currents measured in the two distant SQUID rings. 

In this interdisciplinary work we bridge the gap between electron coherence 
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in mesoscopic physics and nonclassical phenomena in quantum optics. Work on 
entanglement of several mesoscopic devices has been reported in |31j . 

The paper is organized as follows. In section 2 we describe certain one-mode and 
two-mode nonclassical fields, which are relevant to the rest of our work, and derive the 
corresponding Weyl function 32 . In section 3 we discuss the magnetic flux and the 
electromotive force operators, which are the dual quantum variables in our context. 

We subsequently turn our attention to electron interference phenomena. In section 
4 we describe the standard Aharonov-Bohm phase factor in electron interference that 
is induced by a magnetostatic flux. In section 5 we describe the electron phase factor 
operator that is induced by nonclassical electromagnetic fields It is explained that 
the expectation value of the phase factor and, consequently, of the electron intensity 
distribution depend on the quantum state of the external photons. 

We stress that accurate knowledge of the quantum state of the electromagnetic 
field enables us to calculate not only the average intensity of the interfering electrons, 
but also their full statistics (higher order correlations). In section 6 we quantify the 
quantum statistics of the interfering electrons using the autocorrelation function and 
its Fourier transform, the spectral density. It is shown that the quantum statistics 
of the interfering electrons depend on the quantum statistics of the photons |10j . In 
section 7 we describe how two spatially separated electron interference experiments, 
which interact with entangled fields, become correlated 

In section 8 we study the interaction of nonclassical electromagnetic fields with 
mesoscopic SQUID rings. In the case of two distant SQUID rings, which are coupled 
to two entangled electromagnetic fields, we show that the quantum currents tunneling 
through the distant Josephson junctions become entangled [SH]. We conclude with a 
summary of the results in section 9. 

2. Nonclassical electromagnetic fields 

In this section we introduce the nonclassical states of the electromagnetic field that are 
relevant to the rest of our work. We define the Weyl function and provide its value in 
the case of number, coherent, squeezed, and thermal states. It is noted that we use 
theoretical units, in which fc B = h = c = 1. 

2.1. One-mode quantum states of the electromagnetic field 

Nonclassical electromagnetic fields are carefully prepared in a particular quantum state 
and are described by a density matrix p. In this case we know the average values (E), (B) 
of the electric and magnetic fields, the standard deviations AE, AB and also their 
higher moments. Another quantity which describes the fields is the photon counting 
distribution function 



P(N) = (N\p\N). 

Various examples of nonclassical electromagnetic fields are given below. 



(1) 
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2.1.1. Number states The number states \N) are defined as: 



\N) 



|0>. 



Coherent states The coherent states |A) are defined as: 
\A) = D(A)\0) 
where D(A) is the displacement operator 

D(z) = exp(za^ — z*a). 
The photon counting distribution is in this case Poissonian. 



(2) 

(3) 
(4) 



-- exp(— iip)^ 2 + - exp(iLp)a 2 



(5) 



2.1.3. Squeezed states The squeezing operator is defined as 
S(rip) = exp 

where the r, ip are real numbers and r is known as the squeezing parameter. Squeezed 
states | A; rip) are defined by acting on the coherent state | A), with the squeezing operator 

\A;r V ) = S(rip)\A) = S(rip)D(A)\0) . (6) 



In this case P(N) can be sub-Poissonian. The average number of photons is 



(N) sq = 



r 

sinh ( - 
2 
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cosh - 
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1 2 



(7) 



In figure 1 we have plotted the electric field as a function of time in the case of 
coherent and squeezed light. Both the average value (E) and the quantum noise AE are 
shown. The parameters are chosen so that the average value of the electric field is the 
same in both examples. It is seen that the two fields differ in the quantum noise AE. The 
anti-bunching of photons in squeezed states, in comparison to the Poissonian statistics 
in the case of coherent states, is also shown in the figure. These two types of nonclassical 
electromagnetic fields will be used later, in the context of electron interference (i.e., we 
will study the situation where these nonclassical electromagnetic fields are coupled with 
electron interference devices). It will be shown there that they produce different results 
for the electron interference. 

2.1.4- Thermal states The thermal states are defined through the density matrix 
Pth = [1 — exp(— (3uj)] exp(— ftua^d) 

oo 

= [l-exp(-/3w)]£exp(-/3wn)|n)(n| (8) 

?1=0 

where f3 is the inverse temperature. In this case the average number of photons is 

1 
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Figure 1. The electric field of coherent and squeezed light as a function of time. 
Both the average value (E) and the quantum noise AE are shown. The anti-bunching 
of photons in squeezed states, in comparison to the Poissonian statistics in coherent 
states, is also shown. 



2.2. Weyl functions 

The Wigner and Weyl (or characteristic) functions play an important role in quantum 
mechanics |32] • The Weyl function that corresponds to a quantum state described by a 
density operator p is defined in terms of the displacement operator of equation (pEJ) as 

W(z) = Tr[pD(z)\. (10) 

The tilde in the notation reflects the fact that the Weyl function W is the two- 
dimensional Fourier transform of the Wigner function W. The W(z) is a complex 
function, in general, whose absolute value obeys 

0<\W(z)\<l. (11) 

For later use we give the Weyl function for various states. We start with the 
following relation 



(M\D(z)\N) 



NT 
Ml 



1/2 



z M N exp 



■ M—N 



(W 2 ) 



(12) 



where the L% are Laguerre polynomials [SI] . Therefore the Weyl function for a number 
state \N) is 

¥W(M 2 )- (13) 



exp 



The Weyl function for a coherent state \A) is 

I ~|2 



W coh (z) = exp 



+ i2\Az\ sin(arg2; — arg A) 



(14) 
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The Weyl function for a squeezed state \A; rip) is 
W sq (z) =exp(-Y + tX), 

X = 2\Az\ cosh sin(argz — argv4) — sinh \ sin(argz + argv4 + cp) 
\z\ 2 

Y = [cosh(r) + sinh(r) cos(2 arg z + <p)\ . 
Finally for thermal states we have 

W th (z = Ce*"*) = exp 

where ( is a real number. These relations have been given in reference |25| . 



T coth (^ 



(15) 



(16) 



2.3. Two-mode quantum states: separability versus entanglement 

Nonclassical electromagnetic fields with several modes allow for correlations between 
the distinct field modes. The nature of the correlation can be classical or quantum |17j . 

Let p be the density matrix that describes a two-mode nonclassical electromagnetic 
field. Then the density matrices of the two fields are 

p A = Ttb (p) , Pb = Tr A (p). (17) 

The density matrix p for the two-mode electromagnetic field state is factorizable if 
Pfact = Pa <S) Pb- The density matrix p is separable if 

Psep = P kPA,k ® pB,k (18) 



where Pk are probabilities. In all other cases the density matrix p ent is entangled. 



2.3.1. Two-mode number states For later use we consider the (mixed) separable density 
operator 

Pse P = \{\N l N 2 ){N l N 2 \ + \N 2 N l )(N 2 N 1 \). (19) 

We also consider the (pure) entangled state \s) = 2~ 1 / 2 (\NiN 2 ) + \N 2 Ni)), for example. 
The corresponding density operator is 

Pent = Psep + \ {\N X N 2 ) {N 2 N X \ + \N 2 N X )(N X N 2 \). (20) 

Clearly in this example the p sep and the p e nt differ only in the above nondiagonal 
elements. In both the separable and the entangled case the reduced density operators 
of equation (|T7jl are given by 

Psep,A = Pent,A = Pse Pi B = Pent,B = - ( | iVi> <iV"i | + | i\T 2 > <iV 2 1) . (21) 
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2.3.2. Two-mode coherent states We consider the two-mode coherent states in the 
classically correlated state 

p scp = i(|AA 2 )(AA 2 | + lA^X^AI). (22) 

In this case the reduced density operators are 

Pscp,A = Pscp,B = ^(l^iX^il + \A 2 ) (A 2 \). (23) 

We also consider the entangled state \s) = Af(\A 1 A 2 ) + |A 2 Ai)) with density 
operator 

p cnt = 2Af 2 Pscp + M 2 (\A 1 A 2 )(A 2 A 1 \ + I^XA^I) (24) 
where the normalization constant is given by 

M= [2 + 2exp(-|A 1 - A 2 \ 2 )y V2 . (25) 
In this case the reduced density operators are 

Pcnt,A = Pent,B=N 2 (\A 1 )(A 1 \ + \A 2 ){A 2 \+ X \A 1 ){A 2 \+ X *\A 2 ){A 1 \) (26) 



where 



X = (A^A*) = exp (- l -^f - !M + A{A 2 ) . (27) 



3. Magnetic flux operator 

We consider a monochromatic electromagnetic field of frequency iv, at sufficiently low 
temperatures T <C iv, so that the quantum noise is greater than the thermal noise. In 
this case the vector potential A\ and the electric field % are dual quantum variables. 
For a loop C, which is small in comparison to the wavelength of the electromagnetic 
field, the A^Ei are integrated around C and yield the magnetic flux <fi = § c Aidxi and 
the electromotive force Vemf = §c^i^ x %i correspondingly, as dual quantum variables. 
In terms of these variables the photon creation and annihilation operators are 

a f = -^-(0 - ^"Vemf), a = (J> + w^Vemf) , (28) 

where £ is a constant proportional to the area enclosed by C. Consequently the magnetic 
flux operator is 0(0) = 2~ 1 / 2 £(a^ + a) and its evolution in time is given by 

4>(t) = exp(itW)0(O) exp(-im). (29) 

The Hamiltonian H of the system is 

H = uo(a)a + 1/2) + H mt . (30) 

In the external field approximation we ignore the interaction Hamiltonian 7i m t and we 
obtain 

4>{i) = —j= exp(iut)d^ + exp(—iut)d . (31) 
V2 L ' J 
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This is a good approximation when the flux due to back reaction is small in comparison 
to the external flux. 

The expectation value of the flux and the quantum uncertainty A0 are given by 
<0(t)> = Tr[p0(t)], A0 = [<0 2 (t)> - (0(t)) 2 ] 1/2 - (32) 

For example, in the case of number states p = \N)(N\ we get 

1\ 1/2 



(0(t))num = 0, (A0) num = IN + -j . (33) 

For coherent states p = \A)(A\ we have 

(0(t)) coh = 2 1 / 2 |A|cos(^-argA), (A0) coh = 2^ 2 . (34) 
In the case of squeezed states p = \A; rip) {A] rp>\ we obtain 

(0(t))sq = -2 1/2 |A||Z|cos(cjt + argZ), (35) 

(A0) sq = 2^ 1/2 [cosh(r) - sinh(r) cos(2cut + <p)] 1/2 , 

Z = sinh(r/2) exp[z(arg A + (p)] — cosh(r/2) exp(— i arg A). 
Finally for the thermal states of equation (jHj) we get 



(0(t))th = 0, (A#)th = ^ 



cot h ( -/3a; 



1/2 



(36) 



4. Aharonov-Bohm phase factor induced by a magnetostatic flux 

We consider a two-path electron interference experiment, as shown in figure 2. 

C 



1(A) 



Figure 2. Aharonov-Bohm experiment. The electrons follow the lowest winding paths 
Co, C\ in a field-free region. The loop Co — Ci is threaded by a magnetostatic flux $. 

The wavefunctions corresponding to paths Cq and Ci are ipo and ^i, respectively. 
In the presence of magnetic flux $ threading the loop Co — C\, we get the electron 
intensity 

I(x) = |^o| 2 + |^i| 2 + 2|^i| cos(x - e$), (37) 
where x is the phase difference between the two paths: 

x(A) = arg(^o) - arg(^i). (38) 
If we assume equal splitting (i.e., |-?/>o| 2 = 1/2 = |^i| 2 ) then 

I(x) = l + cos(z-e$). (39) 
The visibility of the intensity J, defined as 

(40) 



-^max -^min 



-^max "I - Imin 

is equal to one in this case. 
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5. Phase factor operator induced by nonclassical electromagnetic fields 

In this section we consider a mesoscopic electron interference device (~ 0.1/iin) in a 
microwave waveguide at low temperatures (10 — lOOmK). The electric field is parallel 
to the plane of the electron paths and the magnetic field is perpendicular to it (figure 3). 
The electron intensity is given by equation ()37jl where the flux <3> is now time-dependent. 

We next consider the case where the microwaves are nonclassical. In this case $ is 
a quantum mechanical operator and its expectation value with respect to the density 
matrix p of the microwaves gives the observed electron intensity: 

I(x,t) = l + Tr[pcos(x-e0)] = 1 + $l[e ix W(\)} = 1 + \W(\)\ cos{x - arg[W(A)]}. (41) 

Here W is the Weyl function of the density matrix p defined in equation (fTUj) . and we 
define 

y/2' 

If we compare and contrast equation (}3*9"j) for classical microwaves, with equation (j41|) 
for nonclassical microwaves we see that the visibility is reduced in the second case from 
1 to | W(A) |. This is due to the quantum noise in the nonclassical microwaves as can be 
seen from the expansion 



A = iqexp(iut), q = —=. 



(42) 



l^(A)| 2 = 1 - |[(AX) 2 + (AP) 2 ] - |-[(AX) 2 - (AP) 2 ] cos(2^) - . . . 

where X = £ -1 </>, P = (^0 _1 ^emf, and AX, AP are the corresponding uncertainties. 
The |W"(A)| is less than 1 due to the non-zero values of the quantum noise AX, AP. 

Results are given below for the electron intensity I(x, t) that corresponds to 
irradiation with several quantum states. We choose the point x = 0, for simplicity. 



(43) 



Figure 3. Modified Aharonov-Bohm experiment in the presence of an electromagnetic 
field. The field travels in the waveguide with the magnetic field perpendicular to the 
plane of the electron paths Co, C\ and the electric field parallel to it. 



5.0.3. Number states For the number states of equation (J2J), using equation (|T3j) . we 
get 



(t) = l + exp(-^)lw(g 2 ). (44) 
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Figure 4. Vacuum-induced phase factor for the charges as a function of time, uit, for 
u) = 1CP 4 , corresponding to the case of number states (broken line), coherent states 
(solid line), squeezed states (line of circles), and thermal states (line of stars). The 
average number of photons (N) is zero; the squeezing parameter r is 0.5. In subplot 
(a) the |Vy(A)| is shown and in (b) the arg[W(A)] is shown. 



5.0.4- Coherent states For the coherent states of equation (J3J), using equation ()14jl. we 
get 

I coh (t) = 1 + exp cos \HA cosM - arg A)} . (45) 

In this case the result is very similar to the classical result of equation (}39|) but the 
visibility is slightly reduced from 1 to exp (— ^j. The quantum noise of the coherent 
states slightly destroys the interference and reduces its visibility. Even in the absence of 
microwaves (vacuum state) we get a reduction in the visibility due to the vacuum noise. 

In figure 4 we have plotted the expectation value of the phase factor operator, 
which is induced by the electromagnetic vacuum, the coherent states of equation (JHJ), 
the squeezed states of equation © for r = 0.5, and the thermal states of equation ©. 

5.0.5. Squeezed states For the squeezed states of equation ©, using equation ([To]) , we 
get 

I sq (t) = 1 + exp(-F 1 ) cospfi), (46) 
q 2 

Y\ = — [cosh(r) — sinh(r) cos(2u;t + cp)\, 



X x = 2q\A\ 



cosh (— j cos(cut — arg A) — sinh f-^j cos(tot + arg A + (p) 



We note that in the case of squeezed vacuum (A = 0) the intensity I sq (t) contains all 
the frequencies 2Ku where K is an integer (after a Fourier expansion). In contrast in 
the case of coherent states we get all the frequencies Kuj. The factor of 2 in the case of 
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squeezed vacuum is related with the fact that the squeezed vacuum is a superposition 
of even number states only. Therefore the electrons can only absorb an even number of 
photons (there are no odd number states in this quantum state). In this case the result 
is qualitatively different from the classical result. 

5.0.6. Thermal states For the thermal states of equation (JHJ), using equation (fTflj). we 
have 



J t h(t) = 1 + exp 



(47) 



6. Quantum statistics of the interfering electrons 



There are various quantities that can be used to describe the quantum statistics of the 
interfering electrons. In the previous section we studied the electron intensity and here 
we consider higher order correlations ^Hj. We compare and contrast the results for the 
two cases of classical and nonclassical microwaves. 

6.1. Autocorrelation function of the electron intensity in the case of classical 
microwaves 

In general for a function I(t) the autocorrelation function is defined as 

r(r)= lim J- [ T F(t)I(t + r)dt. (48) 

T^oo II J-T 

The following properties are well known (e.g., see reference |35j): 

r(-r) = r*(r), r(o)>o, |r(r)|<r(o). (49) 

The normalized autocorrelation function is defined as 

An expansion of T(r) into a Fourier series yields the spectral density coefficients 

S K =— T(T)exp(-iKnr)dT (51) 

Z7T JO 

oo 

T(r) = S K exp(iKQr). 

K=~oo 

The property T(— r) = T*(r) of equation (jj^j) guarantees that the coefficients Sk are 
real numbers. If the autocorrelation function is purely real then the spectral density 
coefficients obey the relation Sk = S_k- But if T(r) is complex then, in general, 
Sk 7^ S-k and we refer to this as an asymmetry in the spectral density. 

As an example we consider classical microwaves of frequency uo with magnetic flux 
of the form 

<j>{t) = 0isin(cjt). (52) 
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In this case the electron intensity at the point x = on the screen is 

Ici{t) = 1 + cos[e0i sm(ut)]. (53) 
Therefore the autocorrelation function is 

oo 

r cl (r) = [1 + J o (e0i)] 2 + 2 fc(e</>i)f cos{2Kur), (54) 

K=l 

where the J n (z) are Bessel functions jHlj, and the spectral density coefficients are 

S Q =[1 + J o (e0i)] 2 , S K = [J 2 K(e<pi)] 2 . (55) 




Figure 5. 1(0, t) for the electrons as a function of time, cut, for w = 10 -4 , 
corresponding to irradiation with number states (broken line), coherent states (solid 
line), squeezed states (line of circles), and thermal states (line of stars). We have chosen 
(N) = 17, in all four cases, and r = 4.2 for the squeezed states. For comparison, we 
also show the electron intensity of equation H53|l corresponding to classical microwaves 
(line of crosses), for fa = (2(A^)) 1 / 2 . 



6.2. Autocorrelation function of the electron intensity in the case of nonclassical 
microwaves 

In this case the electron intensity I(t) = 1 + cos[e0(t)] is an operator. Consequently the 
autocorrelation function T(r) of equation (J4"%j) is in this case defined as 

r(r) = lim f T Tr[pP(t)I(t + r)]dt. (56) 

The values of T(r) for the electric charges have been derived for irradiation with 
various nonclassical microwave states in [TO] . m the following we present numerical 
results, which illustrate the electron correlation properties, and allow for a comparison 
between the effects of classical and nonclassical microwaves. 
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Figure 6. Real and imaginary parts of 7(7-) of equation (|5U[) for the electrons as a 
function of time, wr, for uj = 10~ 4 , corresponding to irradiation with number states 
(broken line), coherent states (solid line), squeezed states (line of circles), and thermal 
states (line of stars). The parameters are (N) = 17, r = 4.2, and for the case of 
classical microwaves (line of crosses) we have <pi = ^(iV)) 1 / 2 . 
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Figure 7. Spectral density coefficients Sk of equation l|51(l for the electrons as a 
function of K, corresponding to irradiation with classical microwaves (first column 
from the left), and nonclassical microwaves in number states (second column), coherent 
states (third column), squeezed states (fourth column), and thermal states (fifth 
column). The parameters are (N) — 17, r — 4.2, and tpi = {2{N)) 1 / 2 . 



6.3. Numerical results 

In order to make the comparison meaningful, in the numerical calculations (figures 5-7) 
the number of photons in the number states is equal to the average number of photons 
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in the coherent, squeezed, and thermal states; we have chosen 

N = (AOcoh = (N) sq = (iV>th = 17. (57) 
For comparison with the case of classical microwaves we have chosen the amplitude of 
the classical magnetic flux to be <fii = (2(N)) l l 2 . The frequency of the microwaves is 
uj = 10~ 4 in units where ks — h — c — 1. The squeezing parameter is r = 4.2 and the 
other parameters are £ = 1, &rgA = 0, ip = 0. 

In figure 5 we show the electron intensity I(t) for x = as a function of cut for 
irradiation with number states (broken line), coherent states (solid line), squeezed states 
(line of circles), and thermal states (line of stars). For comparison, we also show the 
electron intensity of equation ()53|) corresponding to classical microwaves (line of crosses). 

In figure 6 we show the real and imaginary parts of j(t) of equation (|5()jl for the 
electrons as a function of ujt corresponding to irradiation with number states (broken 
line), coherent states (solid line), squeezed states (line of circles), and thermal states 
(line of stars); we have also included the results in the case of classical microwaves (line 
of crosses) for comparison. It is seen that different states of the electromagnetic field 
lead to different electron correlation properties. The imaginary part of the normalized 
autocorrelation function vanishes only for irradiation with classical microwaves, but it 
is nonzero for the four cases of nonclassical microwaves. 

In figure 7 we plot the spectral density coefficients Sk of equation (JoTj) for the 
electrons as a function of K, corresponding to irradiation with classical microwaves (first 
column from the left), and nonclassical microwaves in number states (second column), 
coherent states (third column), squeezed states (fourth column), and thermal states 
(fifth column). 

7. Entangled currents in distant electron interference experiments induced 
by entangled photons 

In this section we consider two electron interference devices that are far from each 
other A photon source irradiates the two experiments with correlated two- mode 
nonclassical microwaves. Each microwave field mode is coupled to one of the two 
experiments. The experiment is depicted in figure 8. It will be shown that the photon 
correlations are transferred to the electron interference experiments, which become 
correlated. The nature of their correlation depends on whether the external photons are 
separable (classically correlated) or entangled (quantum mechanically correlated). 

Let p be the density operator describing the two-mode nonclassical electromagnetic 
field. The first mode of frequency iO\ interacts with electrons in experiment A and 
its density matrix is given by p\ = Tre(p). Similarly the second mode of frequency 
uj 2 interacts with electrons in experiment B and its density matrix is pe = T^a{p)- 
The density matrix p can be factorizable (i.e., the field modes are independent of each 
other), separable (the field modes are classically correlated), or entangled (the field 
modes are quantum mechanically correlated). The difference between these cases has 
been discussed in section 2. 
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Figure 8. Two electron interference experiments which are far from each other are 
irradiated with nonclassical electromagnetic fields. The two electromagnetic fields in 
the two experiments are produced by the source Sem and are correlated. 



7.1. Correlations of the electron intensity distributions 

The nonclassical magnetic flux A that influences the electron interference in A gives 
rise to the phase factor operator exp(ie0 A ). This phase factor induces the electron 
intensity distribution / A (£ A ), which is given by 

I A (x A ) = Tr{p A [l + cos(x A - e0 A )]} = 1 + \W(\ A )\ cos{a; A - arg[iy(A A )]} (58) 

where Aa = iq exp(ia>ii). Similarly in experiment B, which is influenced by a nonclassical 
magnetic flux B , one obtains the intensity 

I B (x B ) = Tr{p B [l + cos(x B - e0 B )]} = 1 + \W(\ B )\ cos{x B - arg[W>(A B )]} (59) 

where A B = iqexp(iuj2t). 

The electron intensity on the interference screen of experiment A (or B) is 
calculated by tracing the intensity operator with respect to the density matrix that 
describes the corresponding electromagnetic field mode (p A or p B ). The results, I A (x A ) 
and I B (a; B ), are proportional to the probability of detecting an electron at a point x A in 
A or a point x B in B. It is also possible to define the joint electron intensity I(x A ,x B ), 
which is related to the probability of a simultaneous detection of electrons at x A and 
x B . This joint intensity is controlled by the full density matrix p for the two- mode 
electromagnetic field, that is, 

I(x A , x B ) = Tr{p[l + cos(x A - e0 A )] [1 + cos(x B - e0 B )]}. (60) 

If the two-mode electromagnetic field is factorizable (p = p A <E> p B ) then the 
joint electron intensity J(x A ,x B ) is simply equal to the product / A (x A )J B (x B ) of the 
independent intensities. However, if the two field modes are classically or quantum 
mechanically correlated then this is not true, in general. In order to quantify this we 
can define the ratio 

fi= ,fy°> (en 
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which is equal to one only for independent electron intensities. In other words, whenever 
R takes values not equal to one, the electron intensity in experiment A is correlated to 
the electron intensity in experiment B. In what follows we consider particular examples 
that illustrate the effect. 

7.2. Examples and numerical results 

As an example we consider two-mode separable and entangled microwaves in number 
states[in]. We compare and contrast the effects of the separable state 

p sep = i(|00>(00| + |ll)(ll|) (62) 

and the entangled state 2~ 1//2 ( 1 00) + |11)) with density matrix 

Pent = Psep + ^(|00)(11| + |H)(00|). (63) 

It is noted that the p sep and the p ent differ only in the above nondiagonal elements. The 
reduced density operators that describe the electromagnetic field in A and B are, in 
both cases, 

Psep,A = Pent, A = Psep.B = Pent,B = ~ ( 1 0) (0 1 + 1 1) ( 1 1 ) . (64) 




Figure 9. R scp of equation l|65(l as a function of xa,xb £ [— 2%, 2%]. Here 
min(_R sop ) = 1.0001 and max(_R scp ) = 1.2471. The frequencies are lo\ — 1.2 x 10 -4 and 

The ratio of equation (JBTj) corresponding to the separable state p sep is given by 

„ , s 1 + «(cosxa + COSXr) + 7COSXaCOSX B , nr ^ 

R Se p{X K ,X B ) = — r— r 65 

(1 + acosxAj(l + acosxej 
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where 



a 



2-q 2 



exp -- 



7 



It can easily be shown that 
1 + 2a + 7 



< R sep (x A ,x B ) < 



^exp(-g 2 )[l 



1 - 2a + 7 



(66) 



(67) 



(1 + a) 2 (1-aj 
which in our example leads to min(i? scp ) = 1.0001 and max(i? sep ) = 1.2471. In figure 9 
we plot the R sep (xA, %b) against screen positions xa and xb for microwave frequencies 
u x = 1.2 x 10" 4 and cu 2 = 10" 4 . 




Figure 10. i? on t of equation H68JI as a function of xa, xb G [— 27r, 27r] for (cji+o;2)< = 7r. 
The top and bottom plateaus show the max(i? scp ) = 1.2471 and min(i? sep ) = 1.0001, 
respectively. The frequencies are u>i — 1.2 x 10~ 4 and U02 = 10~ 4 . 



In the case of the entangled state p cnt of equation the ratio is 



2n sin xa sin xb cos[(o; 1 + ^2)*] 



R ent (x A ,XB,t) = R sep (xA,x B ) - q exp(-g )— — — -. (68) 

(1 + a cosxa)(1 + a cosxb) 

The i? ent oscillates in time around the -R sep , with frequency ui\ + and exceeds 

periodically the bounds of the inequality for R sep in ()67jl . In figure 10 we plot the 

-Rent(^A, £b) against screen positions xa and xb for (u\ + Uz)t = it and the same 

microwave frequencies as in the previous figure. The two horizontal surfaces represent 

the min(i? sep ) (bottom plateaux) and the max(i? sep ) (top plateaux). In figure 11 we 

compare the R sep (line of circles) and the R ent (solid line) against time {uj\ + 6^2) i for 

fixed screen positions xa = 0.97T and xb = 1.0257T. 
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Figure 11. Comparison of R en t (solid line) and R scp (line of circles) of equations (1G8I) 
and Ht)5fl . respectively, for xa = 0.9tt and xb = 1.0257T as a function of dimensionless 
time. The frequencies are u)\ = 1.2 x 10~ 4 and ll>2 = 10 -4 . 



8. Interaction of mesoscopic SQUID rings with nonclassical 
electromagnetic fields 

In this section we investigate application of the above ideas in the context of mesoscopic 
superconducting quantum interference device (SQUID) rings. In the first instance we 
introduce mesoscopic SQUID rings and describe how they interact with nonclassical 
electromagnetic fields [23 I2SI- m this case the Josephson currents are quantum 
mechanical operators, whose expectation values with respect to the density operator of 
the external photons, yield the observed currents. Subsequently, we apply the general 
concept described in the previous section to the case of two distant SQUID rings, each 
of which is coupled to a single mode of a two-mode nonclassical electromagnetic field 
|29j . It is shown that the photon correlations are transferred to the Josephson currents 
in the distant superconducting devices. 



8.1. Mesoscopic SQUID ring 

Consider a superconducting ring of mesoscopic area £ < 10~ 8 cm 2 , which is interrupted 
by a Josephson junction (weak link), as shown in the figure below. In this case the 
capacitance C across the Josephson junction is very small and at low temperatures 
T < O.IK the behaviour of the SQUID is nonclassical t 22j in the sense that the Coulomb 
charging energy for a Cooper pair of charge 2e, 

Ec = &f-, (69) 
becomes a significant parameter. 
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Figure 12. Mesoscopic SQUID ring: superconducting ring of mesoscopic dimensions 
is interrupted by a Josephson junction. 

Under these conditions the charge Q = —i(2e)ds through the junction and the phase 
difference across the junction 5 are conjugate operators, which obey the commutation 
relation [5, Q] = i2e. In this case the Josephson current also becomes an operator which 
is a sinusoidal function of the phase difference across the junction, 

I = I c sin 5, (70) 

where J c is the critical current. 

8.2. Interaction with classical microwaves 

We consider a mesoscopic SQUID ring interacting with a monochromatic electromag- 
netic field. The magnetic flux <p(t) is threading the SQUID ring and the phase difference 
across the junction is 5 = 2e<fi(t). 

In the classical case the Josephson current I and the phase difference 5, are classical 
numbers. Therefore for a magnetic flux with a linear and a sinusoidal component, 

(j>{t) = o + Uit + Msin(cJit), (71) 

we get the current 

I = h sin[2e0(t)] = h sin[2e0 o + 2eUit + 2eusin(u; 1 t)]. (72) 

Using the well-known identity 

oo 

exp(ra sin z) = ^ J n (u) exp(inz) , (73) 

n=— oo 

we can easily show that the current can be expanded as 

oo 

I = h Jn(2eu) sin[(2eUi + nw^t + 2&f> ]. (74) 

n=— oo 

Calculating the time-averaged value I^c of the current I we see that when 

2eU = Noo (75) 

where N is an integer, we get 

-fdc = h J-jv(2ew) sin(2e0 o ), (76) 

otherwise the I& c vanishes. These integral values of the voltage are usually referred to 

as Shapiro steps. 
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8.3. Interaction with nonclassical microwaves 

We now study the effect of nonclassical microwaves on the Josephson current of 
a mesoscopic SQUID ring operating at low temperatures. We use the external field 
approximation, and ignore the back-reaction. This is a good approximation when the 
external electromagnetic fields are much stronger than the fields induced by the currents 
circulating the mesoscopic devices. 

We consider the irradiation of a mesoscopic SQUID ring with monochromatic 
nonclassical microwaves of frequency u\. In addition to that the ring is threaded by 
the classical flux O + V\t and the total flux is $(£) = O + V\t + 4>(t). Therefore the 
quantum current is in this case given by 

I A = h sin |2e0 o + 2eV\t + q' exp{iuit)a^ + exp(— iuoit)a j 

= I^s {exp[i(u A t + 2e<po)]D[iq' exp(iuii)]} (77) 

where 

UA = 2eVi, q = V2e£. (78) 

It is noted that the scaled electric charge q' has twice the value of q of equation 
(142)1 . because in this case we have pairs of electrons. The experimentally measured 
current is calculated by tracing with respect to the density operator pa for the external 
electromagnetic fields, that is, 

(I a) = Tt(paIa) = Ii^[exp(iu A t)W((7 A )], cr A = irf exp(iu)it) . (79) 

As an example we consider microwaves in coherent states. For comparison with the 
classical case of equation (|75|l we take coherent states with A = 2~ l / 2 u and argv4 = 0. 
In this case we get Shapiro steps, as in the classical case, but the dc current is now 
reduced by a small factor: 

4c 0h) =exp(-^)/ dc . (80) 

We also consider the case where the microwaves are in a squeezed vacuum. The 
squeezed vacuum is a superposition of even number states only. In this case we get 
even Shapiro steps only. A physical interpretation of this result is that the electrons can 
only absorb an even number of photons (there are no odd number states in this quantum 
state). Similar results can be proved for even Schroedinger cats (Af(\A) + 1 — A)), which 
are superpositions of even number states also. We stress that in this case the result is 
qualitatively different from the classical result, in the sense that the odd Shapiro steps 
are absent. 



8.4- Entanglement of distant mesoscopic SQUID rings 

We consider two spatially separated mesoscopic SQUID rings, which we refer to as A and 
B. They are irradiated with microwaves that are described by a density operator p. The 
microwaves are produced by the same source and are correlated. Photons of frequency 
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u}\ interact with device A;and photons of frequency U2 interact with device B. The 
SQUID rings A and B are also threaded by a classical time-dependent magnetic fluxes 
that increase linearly with time (Vxt and V B t, respectively). The proposed experiment 
is illustrated in figure 13. 





Figure 13. Two distant mesoscopic SQUID rings A and B are irradiated with 
nonclassical microwaves of frequencies oj\ and 0J2, respectively. The microwaves are 
produced by the source Sem and are correlated. Classical magnetic fluxes VpJ, and Vgt 
are also threading the two rings A and B, respectively. 



The observed Josephson current in SQUID ring A or B is given by the expectation 
value of the corresponding current operator, 



(I a) = IiTr(p A sin 5 A ) 
(4) = 7 2 Tr(p B sin 5 B ) 



5 A = 2eU A t + 2e0 A (t), 
5 B = 2eV B t + 2e4> B {t), 



where 



Mt) 



V2 



exp(iuit)di^ + exp(—iuit)di 
exp(iu2t)a2 + exp(— wj%t)a2 



(81) 
(82) 

(83) 
(84) 



in accordance with the formalism developed in section 3 and assuming that both rings 
have the same area £. The (J A ) has been written in terms of the Weyl function PU(cr A ) 
in equation f)79j) : and similarly for B. 

The expectation value of the product of the two current operators is given by: 

(44) = hh^ifi sin S a sin S B ). (85) 

The correlations between the observed electron currents can be quantified by defining 
the ratio 

(44) 



(4) (4)' 



(86) 
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where the superscript (c) indicates that this quantity corresponds to currents. For 
factorizable density matrices pf act = Pa ® Pb we easily obtain the ratio R^. t = 1, 
identically. For separable density matrices p sep of equation (|T%j) we get 

i?S = — gigiiMM, . (87) 

P (E*P*</a*»(EiPi</bi» 
We also calculate the higher moments of the currents 

(J A ) =/ 1 2 Tr[p A (sina A ) 2 ], (88) 

(/|> =/ 2 2 Tr[p B (sin ( TB) 2 ], (89) 

(IVb) = A% 2 Tr[p(sina A ) 2 (sina B ) 2 ]. (90) 

The expectation value (Ia^b) quantifies the quantum statistics of the electron pairs 
tunneling the junctions in the two SQUID rings. Consequently the ratio 

R (c2) = (i 2 Jl) (91) 

is a measure of the photon-induced correlations of the quantum statistics of the tunneling 

(c2) 

electrons. For factorizable density matrices we easily see that -Rf act = 1- For separable 
density matrices we get 

^>(c2) _ ^iPi(^Ai)(^Bi) ^g2) 



sep 



8.5. Examples and numerical results 

We present examples in which we compare and contrast the influence of a classically 
correlated two-mode microwave state with a quantum mechanically correlated one, on 
the Josephson currents. The two-mode microwaves are in both number and coherent 
states. 

8.5.1. Number states Firstly we consider the separable density operator p sep of equation 
(ITTjj) and the entangled density operator p cnt of equation (|271jl for number states. 
For the p sep of equation (fT9"j) we calculate the currents in A and B: 

(I A ) = hC sm(uj A t), (93) 
(J B ) = hC sin(u B t), (94) 

Co = ^exp (-*Pj [L Nl (q' 2 ) + L N2 (q'% (95) 

where the L"(x) are Laguerre polynomials [33]. It is noted that in this case the currents 
are independent of the microwave frequencies ui,U2- The second moments 
of the currents in A and B, defined by equations (fHHJ) and (jEHJ), respectively, have also 
been calculated: 

(^) = y[ 1 -Cicos(2^ A t)], (96) 
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(Ib) = S[i-Cicos(2wBf)], 
Ci = iexp(-2g' 2 )[L^(4g' 2 ) + L W2 (4g' 2 )]. 
The expectation value of the product of the two currents is 

(lA^sep = hhC 2 sin(a; A t) sin(o; B i), 
C 2 = exp(-q' 2 )L Nl (q' 2 )L N2 (q' 2 ). 

Consequently the ratio of equation (jHBj) is 
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(97) 
(98) 

(99) 
(100) 

Co 2 ~ [LnM 2 ) + L N2 (q' 2 )r (101) 

In this example the is time- independent; it depends only on the number of photons 
Ni, N 2 in the two-mode microwave field. 



R(c) = C, 4L Nl (q' 2 )L N2 (q 2 ) 
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Figure 14. Comparison of i?sop for the separable number state of equation (|19f) (line 
of circles) and Ri% for the separable coherent state of equation 11' I'll (solid line) for 
N% = 1,N 2 = 3 and Ai = 1,A 2 = 3 1/2 as a function of {ujx—uj^jt, where Wi = 1.2xl0~ 4 
and ui2 = 10 -4 . 

For the p ent of equation the (/a), (^b) are the same with those presented in 
equations (jHSJ), (jHU); and the (Ij}, (Iq) are the same as in equations (JHBJ), (}9Tj) . However 
the {I\Ib) is in this case 

(^A-^B)ent = (^A^sep + -Across ; 

(102) 

where 

/cross = - /l/ 2 C 3 [cOs(cJAt + 0U B t) - (-l) N ^ Na COs(uJ A t - U B t)} COs(Ot), (103) 

^3 =iexp(-g' 2 )LS^(g' 2 )^r Af2 (^)- (104) 
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The /cross causes the (lA-^B)ent to oscillate in time around the (/A-^sep with frequency 

Q = {N 1 -N 2 ){u; 1 -u; 7 ). (105) 

We note that the term J cr0 ss is induced by the nondiagonal elements of p cnt of equation 
P0|). and depends on the photon frequencies ui,u 2 . This term quantifies the difference 
between the effect of separable and entangled microwaves on the Josephson currents. 



(a) NUMBER, N 1= 1, N 2 =3 




DC 

I 

Q- 
<D 

DC 



-0.2 Q 

1 



-1 



(b) COHERENT, A,=1, A 2 =3 



4 

1/2 



(a-a> 2 )X 



Figure 15. The difference iisep— -Rent corresponding to (a) the separable and entangled 
number states of equations (|19|l . I|2U|I : and (b) the separable and entangled coherent 
states of equations (J22J), ijUJl, for Ni = 1, N2 = 3 and Aj_ = 1,A 2 = 3 1 / 2 as a function 
of (u>i — u)2)t, where lj% = 1.2 x and U2 — 10~ 4 . 



In the entangled case the ratio R( c > of equation ((^Bj) can be simplified in two distinct 
expressions according to whether the difference Ni — N 2 is even or odd. In the case 
Ni — N 2 = 2k, the ratio is 



R 



(c) 

ent,2fc 



4L^(g' 2 )L^ 2 (g' 2 ) 



R {c) + 



N 2 



COS (fit). 



(106) 



It is seen that the -R e nt,2fc oscillates around the R^ with frequency f2 given by equation 
(|105|) . If there is no detuning between the nonclassical electromagnetic fields, i.e. 
Qnt2fc 1S constant, although it is still R^t / ■•^ v - 
2k + 1 the ratio is 



iVi 



= u 2 , then i?i 

N 2 = 



R 



(c) 

ent,2fc+l 



(c) 

In both cases the R eT [ t is time- dependent and it is a function of the photon frequencies 
u>i,u> 2 , in contrast to the case of (which is time-independent). 



Ri c l 



7^ -Rsep- m the C£LSe 

cos(fit) 



sop 



[L Nl (q' 2 ) + L^g' 2 )] 2 tan(o; A t) tan(a; B t) ' 



(107) 
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Figure 16. Difference of currents defined by equations (|79|l and (|88|) that are 
induced by separable and entangled photons in coherent states (for SQUID ring A), (a) 
(^A)sep — (^A)ent and (b) (i/^scp ~ (^1)™* corresponding to irradiation with separable 
and entangled coherent states of equations (|22|l and 124|) . for A\ = 1, A% = 3 1 / 2 as a 
function of (uii — uj2)t, where lo\ — 1.2 x 10~ 4 and 0J2 = 1CP 4 . 



(a) NUMBER, N =1, N =3 
0.04 1 1 1 — i — 




Figure 17. Difference of the product of currents (/a^b)sc P — (^A^B)cnt, defined by 
equation (JSSJ, that are induced by (a) the separable and entangled number states of 
equations i|19|). (J2UJ); and (b) the separable and entangled coherent states of equations 
(|2lJ) . for Ni = 1,N 2 = 3 and Ai = 1,A 2 = 3 1/2 , respectively, as a function of 
(wi — W2)i, where = 1.2 x 10~ 4 and ui2 = 10~ 4 . 
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8.5.2. Coherent states We consider the separable density operator p sop of equation (J22J) 
and the entangled density operator p cnt of equation (J2IJ). 

For the separable state of equation (1221) the currents in A and B are 



(/a)so P = y ex P( — ^-){sin[^ A i + 2g'|Ai| cos(^it - 9 X )\ 
+ sin[cj A t + 2q'\A 2 \ cos(to>it - 6 2 )]}, 



108) 



J 2 g'2 

(^b)scp = y ex P( — ^-){sin[cu B ^ + 2g cos(u; 2 t - #i)] 

+ sm[u B t + 2q'\A 2 \cos(u 2 t-9 2 )}}, (109) 

where 6*i = a.rg(Ai), and # 2 = arg(l 2 ). The expectation values of the product of the 
currents, and hence the ratios and Rj$, have been calculated numerically 
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Figure 18. The difference i?sep' — -fi^n? corresponding to (a) the separable and 
entangled number states of equations 1)19(1. I(2Q|I : and (b) the separable and entangled 
coherent states of equations (EH, for iVi = 1, iV 2 = 3 and A\ = 1, = 3 1 / 2 as a 
function of — 0^2 where a>i = 1.2 x 10~ 4 and lo 2 = 1CP 4 . 



For the entangled state of equation ()24|) the current in A is 
(J A ) cnt = 2A/" 2 (/a)scp + Af 2 EF 1 exp f-^) I l5 



(110) 



where 



and 



E = exp[-|li| 2 - \A 2 \ 2 + 2\A X A 2 \ cos(0i - 2 )], 



Fi = {exptglAl^^t) -g|A 2 |S A2 (t)] +exp[-g / |Ail^,i(t) + gV 2 |5 A , 2 (t)]} 
x sinfc^t + ^lAilC^iCt) + g'|l 2 |Cl i2 (t)]. 



(Ill) 



(112) 
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The terms entering the factor Fx are trigonometric functions of the form 
S Ajl = sin(cjit - 8x), S A ,2 = sin(cuit - 6 2 ), 

Ca,i — cos(uxt — 9x), Ca,2 = cos(u;it — 2 ). (113) 
Similarly the current in SQUID ring B is 

(J B ) c „t = 2A/" 2 (/ B >scp + ■N" 2 EF 2 exp (-^j I 2 , (114) 

where 

F 2 = {exv[q\Ax\S B ,i(t) - q\A 2 \S B , 2 (t)} + exp[-q'\Ax\S B ,x (t) + q'\A 2 \S B , 2 (t)]} 

x sm[co B t + q'\Ax\C B> x(t) + q'\A 2 \C B , 2 (t)}, (115) 

and 

S B ,i = sin(cu 2 t - Ox), S B ,2 = sin(uj 2 t - 2 ), 

C Bjl = cos((jj 2 t — 9x), C B ,2 = cos(u; 2 t - 2 ). (116) 

The expectation values of the product of the currents, and hence the ratios R^ t and 
i?ent , have been calculated numerically. 

8.5.3. Numerical results In figures 14-18 we plot the results against dimensionless time 
(ojx — ou 2 )t, where the photon frequencies are lox = 1.2 x 10~ 4 and uo 2 = 10~ 4 . Other 
fixed parameters are the number of photons in the number states: iVi = 1, N 2 = 3; and 
the average number of photons in the coherent states: Ax — 1, A 2 — 3 1 / 2 (we take these 
values so that the microwaves in number and coherent states contain the same average 
number of photons). 

In figure 14 we present the R^ p for the separable number state of equation (fT^j) (line 
of circles) and the for the separable coherent state of equation (}22|) (solid line). It 
is seen that separable photons in different quantum states induce different correlations 
R^ between the Josephson currents in the distant SQUID rings. 

In figure 15 we show the difference R$ — R^t corresponding to (a) the separable 
and entangled number states of equations (JIHJ), P0|): and (b) the separable and entangled 
coherent states of equations (j22j) , l n this case the separable and entangled photons 
induce different correlations R^ between the Josephson currents. 

In figure 16 we present (a) (/ A )sc P - (-Ment and (b) (/ A ) sep - (^i)ent corresponding 
to irradiation with separable and entangled coherent states of equations (|2*2*|) and (J21J). 

In figure 17 we show the (/A-^B)se P — (^A-^B)ent that are induced by (a) the separable 
and entangled number states of equations (flH|) . (f2T)j) : and (b) the separable and entangled 
coherent states of equations (|2^j) . (jUJ). 

In figure 18 we plot the difference R^ — R^t corresponding to (a) the separable 
and entangled number states of equations (fTT?|) . (J211J) : and (b) the separable and entangled 
coherent states of equations ([2^1. (J2IJ). 
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9. Discussion 

We have studied electron interference in mesoscopic devices in the presence of 
nonclassical electromagnetic fields. The phase factor is in this case a quantum 
mechanical operator, whose expectation value with respect to the density matrix of the 
electromagnetic field determines the electron interference. We have presented various 
examples, which show that the quantum noise of the photons destroys slightly the 
electron interference fringes. Related is also the fact that the photon statistics affects 
the interfering electrons. These ideas have also been applied in the context of mesoscopic 
SQUID rings. 

In certain cases we get novel quantum phenomena with no classical analogue. 
For example, in the case of a mesoscopic SQUID ring irradiated with microwaves in 
a squeezed vacuum state we get Shapiro steps only at even multiples of the fundamental 
frequency. 

An important feature of nonclassical electromagnetic fields is entanglement. We 
have considered two distant mesoscopic electron interference devices that are irradiated 
with a two-mode nonclassical electromagnetic field. Each field mode is coupled to one of 
the mesoscopic devices. For entangled electromagnetic fields, the electric currents and 
their higher moments become correlated. 

All our results have been derived within the external field approximation where 
the back reaction (additional flux created by the electrons) is negligible. This is a valid 
approximation in devices with small inductance. 

Most of the experimental work on mesoscopic devices has studied their interaction 
with classical electromagnetic fields, until recently |24j . Our results show that there 
is merit in having a full quantum system where both the mesoscopic device and the 
electromagnetic field are quantum mechanical. In this case we can have purely quantum 
phenomena, without classical analogue, such as the entanglement of distant mesoscopic 
devices. 
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